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Abstract. We give an introduction crystalline color superconductivity, arguing that it is likely to 
' occur wherever quark matter in which color-flavor locking does not occur is found. We survey the 

O . properties of this form of quark matter, and argue that its presence in a compact star may result in 

pulsar glitches, and thus in observable consequences. However, elucidation of this proposal requires 
\ an understanding of the crystal structure, which is not yet in hand. 



INTRODUCTION 

> 

At asymptotic densities, the ground state of QCD with three quarks with equal masses 
is expected to be the color-flavor locked (CFL) phase [1, 2, 3]. This phase features a 
condensate of Cooper pairs of quarks which includes ud, us, and ds pairs. Quarks of all 
colors and all flavors participate in the pairing, and all excitations with quark quantum 
numbers are gapped. 

The CFL phase persists for unequal quark masses, so long as the differences are 



c3 



not too large [4, 5]. It is very likely the ground state for real QCD, assumed to be in 
equilibrium with respect to the weak interactions, over a substantial range of densities. 
(D \ In this phase, chiral symmetry is broken via the locking of left-flavor and right-flavor 

symmetries to color [1]. Terms of order in the effective Lagrangian for the resulting 
pseudo-Goldstone bosons [6] may rotate the CFL condensate in the ^° direction [7], 
resulting in further pseudo-Goldstone bosons [8]. Throughout the range of parameters 
over which the CFL phase exists as a bulk (and therefore electrically neutral) phase, it 
consists of equal numbers of u, d and s quarks and is therefore electrically neutral in 
the absence of any electrons [9]. The equality of the three quark number densities is 
enforced in the CFL phase by the fact that this equality maximizes the pairing energy 
associated with the formation of ud, us, and ds Cooper pairs. This equality is enforced 
even though the strange quark, with mass m s , is heavier than the light quarks. If higher 
order effects do in fact introduce an additional AT condensate, the conclusion that the 
CFL phase is electrically neutral in the absence of electrons remains, because K° mesons 
are neutral. 

If one imagines increasing m s (or, more physically, decreasing y) color-flavor locking 
is maintained until a transition to a state of quark matter in which some quarks remain 
ungapped. This "unlocking transition", which must be first order [4, 5], occurs when 
nr s ~ 4/jAo [9, 10]. In this expression, Ao is the BCS pairing gap, estimated in both 
models and asymptotic analyses to be of order tens to 100 MeV [3]. The strange 
quark mass m s is a density-dependent effective mass [11]. For /u ~ 400 — 500 MeV, 



corresponding to quark matter at densities which may arise at the center of compact 
stars, m s is certainly significantly larger than the current quark mass, and its value is 
not known. In drawing the QCD phase diagram, therefore, there are two possibilities. 
As a function of decreasing fi, one possibility is a first order phase transition directly 
from color- flavor locked quark matter to hadronic matter, as explored in Ref. [10]. The 
second possibility is an unlocking transition [4, 5] to quark matter in which not all quarks 
participate in the dominant pairing, followed only at a lower n by a transition to hadronic 
matter. We assume the second possibility here, and explore its consequences. 

In quark matter in which CFL pairing involving all quarks does not occur, it is likely 
that up and down quarks continue to pair. In this 2SC phase, which was the earliest 
color superconducting phase to be studied [12], the attractive channel involves the 
formation of Cooper pairs which are antisymmetric in both color and flavor, yielding 

a condensate with color (greek indices) and flavor (latin indices) structure {q^qfy ~ 
£ a ^£ a P 3 . This condensate leaves five quarks unpaired: up and down quarks of the third 
color, and strange quarks of all three colors. Because the BCS pairing scheme leaves 
ungapped quarks with differing Fermi momenta, it is natural to ask whether there 
is some generalization of the pairing ansatz, beyond BCS, in which pairing between 
two species of quarks persists even once their Fermi momenta differ. Crystalline color 
superconductivity is the answer to this question. 

THE CRYSTALLINE COLOR SUPERCONDUCTING STATE 

To date, crystalline color superconductivity has only been studied in the simplified model 
context in which one considers only pairing between massless up and down quarks 
whose Fermi momenta we attempt to push apart by turning on a chemical potential 
difference [15, 16, 17, 18], rather than considering CFL pairing in the presence of quark 
mass differences. That is, we introduce 

/u u = /J-b/j 

/j d = £/ + 6>. (1) 

If d/j is nonzero but less than some 8/ii, the ground state is precisely that obtained for 
8/j = [19, 20, 15]. In this 2SC state, red and green up and down quarks pair, yielding 
four quasiparticles with superconducting gap Ao [12]. Furthermore, the number density 
of red and green up quarks is the same as that of red and green down quarks. As long as 
d/u is not too large, this BCS state remains unchanged (and favored) because maintaining 
equal number densities, and thus coincident Fermi surfaces, maximizes the pairing and 
hence the gain in interaction energy. As 8/j is increased, the BCS state remains the 
ground state of the system only as long as its negative interaction energy offsets the 
large positive free energy cost associated with forcing the Fermi seas to deviate from 
their normal state distributions. In the weak coupling limit, in which Aq//j<^ 1, the BCS 
state persists for 8/j < 8/ji = Ao/a/2 [19, 15]. These conclusions are the same (as long 
as Ao/fi <C 1) whether the interaction between quarks is modeled as a point-like four- 
fermion interaction or is approximated by single-gluon exchange [3]. The loss of BCS 
pairing at 8/j = 8/ji is the analogue in this toy model of the unlocking transition. 



If is too large, no pairing between species is possible. The transition between 
the BCS and unpaired states as 8/j increases past dm has been studied in electron 
superconductors [19] and QCD superconductors [4, 5, 20] assuming that no other state 
intervenes. However, there is good reason to think that another state can occur. This is the 
"LOFF" state, first explored by Larkin and Ovchinnikov[13] and Fulde and Ferrell[14] 
in the context of electron superconductivity in the presence of magnetic impurities. They 
found that near the unpairing transition, in a range <8fi< 8fX2, it is favorable to form 
a state in which the Cooper pairs have nonzero momentum. This generalization of the 
pairing ansatz (beyond BCS ansatze in which only quarks with momenta which add to 
zero pair) is favored because it gives rise to a region of phase space where each of the 
two quarks in a pair can be close to its Fermi surface, and such pairs can be created at 
low cost in free energy. Condensates of this sort spontaneously break translational and 
rotational invariance, leading to gaps which vary periodically in a crystalline pattern. 
If in some shell within the quark matter core of a neutron star (or within a strange 
quark star) the quark number densities are such that crystalline color superconductivity 
arises, rotational vortices may be pinned in this shell, making it a locus for glitch 
phenomena [15]. 

In Ref. [15], we have evaluated the width of the window < 8/j < 8^2 for which 
crystalline color superconductivity occurs, upon making the simplifying assumption that 
quarks interact via a four-fermion interaction with the quantum numbers of single gluon 
exchange. 

In the LOFF state, each Cooper pair carries momentum 2q with |q| rs 1.28/j. The 
condensate and gap parameter vary in space with wavelength 7i/|q|. Although the mag- 
nitude |q| is determined energetically, as we sketch below, the direction q is chosen 
spontaneously. The LOFF state is characterized by a gap parameter A and a diquark 
condensate, but not by an energy gap in the dispersion relation: we obtain the quasipar- 
ticle dispersion relations[15] and find that they vary with the direction of the momentum, 
yielding gaps that vary from zero up to a maximum of A. The condensate is dominated 
by the regions in momentum space in which a quark pair with total momentum 2q has 
both members of the pair within ~ A of their respective Fermi surfaces. These regions 
form circular bands on the two Fermi surfaces. Choosing a single, fixed, q means that 
only one circular band on each Fermi surface participates in the pairing. In all work 
published to date, we assume that all Cooper pairs make the same choice of direction 
q. Making this ansatz corresponds to choosing a single circular band on each Fermi 
surface. In position space, it corresponds to a condensate which varies in space like 

(^(x)^(x))ocAe 2 ^ x . (2) 

This ansatz is certainly not the best choice. If a single plane wave is favored, why not 
two? That is, if one choice of q is favored, why not add a second q, with the same |q| but a 
different q? If two are favored, why not three? This question, namely the determination 
of the favored crystal structure of the crystalline color superconductor phase, remains 
open. Note, however, that if we find a region 8/ji < 8/j < 8^2 in which the simple LOFF 
ansatz, with a single q, is favored over the BCS state and over no pairing, then the 
LOFF state with whatever crystal structure turns out to be optimal must be favored in 
at least this region. Note also that even the single q ansatz, which we use henceforth, 



breaks translational and rotational invariance spontaneously. The resulting phonon has 
been analyzed in considerable detail in Ref. [17]. It will be very interesting to use these 
methods to analyze the phonons in more complicated crystal structures. 

Having simplified the interaction by making it pointlike, and simplifed the ansatz by 
assuming the condensate varies like a plane wave, in Ref. [15] we give the ansatz for the 
LOFF wave function, and by variation obtain a gap equation which allows us to solve 
for the gap parameter A, the free energy and the values of the diquark condensates which 
characterize the LOFF state at a given 8/j and |q|. We then vary |q|, to find the preferred 
(lowest free energy) LOFF state at a given d/u, and compare the free energy of the LOFF 
state to that of the BCS state with which it competes. 1 

Crystalline color superconductivity is favored for 5/ji < d/j < 8,1/2- As 8/j increases, 
one finds a first order phase transition from the ordinary BCS phase to the crystalline 
color superconducting phase at d/j = dfi\ and then a second order phase transition at 
8/j = 8/J2 at which A decreases to zero. Because the condensation energy in the LOFF 
phase is much smaller than that of the BCS condensate at 8/j = 0, the value of 8/ji is 
almost identical to that at which the naive unpairing transition from the BCS state to 
the state with no pairing would occur if one ignored the possibility of a LOFF phase, 
namely 8/ji = Ao/\/2. For all practical purposes, therefore, the LOFF gap equation is 
not required in order to determine 8/ji . The LOFF gap equation is used to determine 8/J2 
and the properties of the crystalline color superconducting phase [15]. 

We find that the LOFF gap parameter decreases from 0.23Ao at d/u = d/ui to zero 
at 8/j = b/J2 [15]. The critical temperature above which the LOFF state melts is T c = 
A^/3/2% 2 [21, 16], which means that, except for very close to 8^2, T c will be much 
higher than typical neutron star temperatures. 

In the limit of a weak four-fermion interaction, G — > 0, the crystalline color supercon- 
ductivity window is bounded by 8/ji = Aq/\/2 and 8^2 = 0.754Ao, as first demonstrated 
in Refs. [13, 14]. These results have been extended beyond the G — > limit in Ref. 
[15]. Note that the BCS gap Ao increases monotonically with G. We may therefore use 
Ao to parametrize the strength of the interaction G. This proves convenient because, as 
we have seen, both 8/ji and 8/J2 are given simply in terms of the physical quantity Ao. 
(Writing them in terms of the model-dependent parameters G and A requires unwieldy 
expressions.) As first recognized by Refs. [13, 14], at any fixed 8/j the LOFF phase only 
occurs when the interaction strength (that is, Ao) lies in a specified window. LOFF pair- 
ing does not survive the weak-coupling (Ao — > 0) limit at fixed 8/j, because in this limit 
the width of the band on the Fermi surface in which pairing occurs goes to zero. On the 
other hand, if one takes 8/j and Ao both to zero while keeping d/u/Ao fixed and in the 
appropriate range, LOFF pairing persists down to arbitrarily weak coupling [13, 14, 15]. 



Our model Hamiltonian has two parameters, the four-fermion coupling G and a cutoff A. We often 
use the value of Ao, the BCS gap obtained at 8/j — 0, to describe the strength of the interaction: small Ao 
corresponds to small G. When we wish to study the dependence on the cutoff, we vary A while at the same 
time varying the coupling G such that Ao is kept fixed. We find that the relation between other physical 
quantities and Ao is reasonably insensitive to variation of A. 



OPENING THE CRYSTALLINE COLOR SUPERCONDUCTIVITY 

WINDOW 



The first part of this section is more technical than the rest of this paper. Some readers 
may wish to skip to the text after Eqs. (15,16). 

The variational derivation of the gap equation for the crystalline color superconduct- 
ing phase is somewhat cumbersome [15]. One constructs a variational ansatz in which 
only quarks within a "pairing region" are allowed to pair, minimizes the free energy 
with respect to all variational parameters (two per mode in momentum, color, flavor and 
spin space), and obtains a self-consistency relation which may then be solved to obtain 
A. The intricacy arises from the fact that the definition of the boundary of the pairing 
region involves A itself. In Ref. [16], we provide a diagrammatic rederivation in which 
one simply makes an ansatz for the quantum numbers of the condensate and then "turns 
a field-theoretical crank" and sees this intricate result emerge. We then use the diagram- 
matic derivation to analyze the LOFF phase at nonzero temperature, obtaining the result 
for T c given in the previous section. The diagrammatic formalism also allows us to go 
beyond a point- like interaction, and treat the exchange of a propagating gluon [18]. Be- 
fore presenting the results of this generalization of the interaction, we first sketch the 
diagrammatic derivation of the gap equation. 

In the crystalline color superconducting phase, the condensate contains pairs of u and 
d quarks with momenta such that the total momentum of each Cooper pair is given by 
2q, with the direction of q chosen spontaneously. As noted above, wherever there is an 
instability towards (2), we expect the true ground state to be a crystalline condensate 
which varies in space like a sum of several such plane waves with the same |q|. 

In order to describe pairing between u quarks with momentum p + q and d quarks 
with momentum — p + q, we must use a modified Nambu-Gorkov spinor defined as 



Note that by q we mean the four- vector (0,q). The Cooper pairs have nonzero total 
momentum, and the ground state condensate (2) is static. The momentum dependence 
of (3) is motivated by the fact that in the presence of a crystalline color superconducting 
condensate, anomalous propagation does not only mean picking up or losing two quarks 
from the condensate. It also means picking up or losing momentum 2q. The basis 
(3) has been chosen so that the inverse fermion propagator in the crystalline color 
superconducting phase is diagonal in p-space and is given by 
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where A = YoA^yo and, here, A is a matrix proportional to Cys£ a $ 3 . Note that the 
condensate is explicitly antisymmetric in flavor. 2p is the relative momentum of the 
quarks in a given pair and is different for different pairs. In the gap equation below, we 
shall integrate over po and p. As desired, the off-diagonal blocks describe anomalous 
propagation in the presence of a condensate of diquarks with momentum 2q. 

We obtain the gap equation by solving the one-loop Schwinger-Dyson equation given 

by 

s-\k, q )-s \k,q) = - 8 2 f ^r;% g )rX v B (^-p), 



(5) 



(6) 



where D^ B = D mV 8ab is the gluon propagator, S is the full quark propagator, whose 
inverse is given by (4), and So is the fermion propagator in the absence of interaction, 
given by S with A = 0. The vertices are defined as follows: 

/ y p X A /2 \ 

YaA A / 2 

-(Y^ A /2) r 

oo o -(y,i A /2) T 

In Refs. [15, 16], we introduce a point-like interaction by replacing g 2 D fjV by g*™ times 
a constant G. After some algebra (essentially the determination of S given specified 
above), and upon suitable projection, the Schwinger-Dyson equation (5) reduces to a 
gap equation for the gap parameter A given (in Euclidean space) by 
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where w = |p| 2 — |q| 2 — (ipo + bp) 2 + p 2 + A 2 . We show in Ref. [16] that upon neglecting 
the (small) contributions of antiparticle pairing, this gap equation simplifies to 
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(2tt) 4 (po-^i(p))(po + ^2(p)) 
where £1,2 (p) are defined as in Ref. [15]: 

£i(p)= +6 A /+i(|p + q|-|p-q|) + iy / (|p + q| + |p-q|-2 A /) 2 + 4A 2 sin 2 | 
£ 2 (p)= -5 A /-i(|p + q|-|p-q|) + i V / (|p + q| + |p-q|-2 A /) 2 + 4A 2 sin 2 | (9) 

and |3 is defined as the angle between the up quark momentum q + p and the down quark 
momentum q — p. Upon doing the po integration, we obtain the gap equation derived 
variationally in Ref. [15]: 
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where the "pairing region" T in p-space is given by 



T = {p |#i(p) >0and£ 2 (p) >0} . (11) 

Thus, an exercise in residue calculus has reproduced the blocking regions, excluding 
from the gap equation those regions in momentum space where £i(p) or ^(p) is 
negative. Note that because £i(p) + 2i2(p) > 0, as can be seen from the definitions (9), 
there is no value of p for which both E\ and £2 are negative. Note also that the gap 
equation is dominated by those regions in momentum space where E\ (p) + £2(p) is as 
small as possible, where the integrand in (10) is of order 1/A. These values of p are 
such that both members of a LOFF pair have momenta close to (within ~ A of) their 
respective Fermi surfaces. That is, |p + q| is within A of \i u and |— p + q| is within A 
of [id. The results described in the previous section all follow from analysis of the gap 
equation (10) [13, 14,21, 15]. 

In Ref. [18], we begin with the Schwinger-Dyson equation (5) but this time keep 
the gluon propagator. That is, we analyze the crystalline color superconducting phase 
upon assuming that quarks interact by the exchange of a medium-modified gluon, as 
is quantitatively valid at asymptotically high densities. The medium-modified gluon 
propagator is given by 

D ^7%> + ¥%)-^< <12) 

where £ is the gauge parameter, G(p) and F(p) are functions of po and |p|, and the 
projectors are defined as follows: 

Pfj = 5 (7 -p t pj, P r = P& = 0, P^ = -g, v + Ef^L -Pj v . (13) 

The functions F and G describe the effects of the medium on the gluon propagator. If we 
neglect the Meissner effect (that is, if we neglect the modification of F(p) and G(p) due 
to the gap A in the fermion propagator) then F(p) describes Thomas-Fermi screening 
and G{p) describes Landau damping and they are given in the hard dense loop (HDL) 
approximation by [22] 

F(,) = ,„^(l-^eo(^)), wUhCb W = I|ogg±^ 
\ ( / • \ 2 \ / • \ • 

(14) 



G{P) = / T-T 
^ IPI 



where m 2 = g 2 fi 2 / K 2 is the Debye mass for Nf = 2. The further modification to the gluon 
propagator due to the Meissner effect in spatially uniform color superconducting phases 
has been the subject of much recent work [23], but the Meissner effect in the crystalline 
color superconducting phase has not yet been analyzed. Fortunately, in our calculation 
of d/J2 we we shall only need to study the crystalline color superconducting phase in the 
limit in which A — > 0, and in this limit the expressions (12) and (14) are valid. 



Upon neglecting the (small) contributions from antiparticle pairing, the gap equation 
becomes 



A(*o) = 
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with £i and £2 as in (9), and where cosoc(fc,p) = (k — q) ■ (k — p) and cos|3(fc,p) = 

(q + k) • (q — p). In Ref. [18], we use this gap equation to obtain 8,1/2, the upper boundary 
of the crystalline color superconductivity window. This analysis is controlled at asymp- 
totically high densities where the coupling g is weak. 

At weak coupling, quark-quark scattering by single-gluon exchange is dominated by 
forward scattering. In most scatterings, the angular positions of the quarks on their re- 
spective Fermi surfaces do not change much. As a consequence, the weaker the coupling 
the more the physics can be thought of as a sum of many (1 + 1) -dimensional theories, 
with only rare large-angle scatterings able to connect one direction in momentum space 
with others [24]. In the LOFF state, small-angle scattering is advantageous because it 
cannot scatter a pair of quarks out of the region of momentum space in which both 
members of the pair are in their respective rings, where pairing is favored. This means 
that it is natural to expect that a forward- scattering-dominated interaction like single- 
gluon exchange is much more favorable for crystalline color superconductivity than a 
point- like interaction, which yields s- wave scattering. 

Suppose for a moment that we were analyzing a truly (1 + 1) -dimensional theory. 
The momentum- space geometry of the LOFF state in one spatial dimension is qualita- 
tively different from that in three. Instead of Fermi surfaces, we would have only "Fermi 
points" at ±/u u and ±/^. The only choice of |q| which allows pairing between u and d 



quarks at their respective Fermi points is |q| = 8/j. In (3 + 1) dimensions, in contrast, 
|q| > d/j is favored because it allows LOFF pairing in ring-shaped regions of the Fermi 
surface, rather than just at antipodal points [13, 14, 15]. Also, in striking contrast to the 
(3 + 1) -dimensional case, it has long been known that in a true (1 + 1) -dimensional the- 
ory with a point-like interaction between fermions, 8^2/Ao — > 00 in the weak-interaction 
limit [25]. 

We expect that in (3 + 1) -dimensional QCD with the interaction given by single- 
gluon exchange, as /u — > °° and g(/S) ^ the (1 + 1) -dimensional results should be 
approached: the energetically favored value of |q| should become closer and closer to 8/j, 
and 8//2/A0 should diverge. We derive both these effects in Ref. [18] and furthermore 
show that both are clearly in evidence already at the rather large coupling g = 3.43, 
corresponding to fi = 400 MeV using the conventions of Refs. [26, 27]. At this coupling, 
8^2/Ao ~ 1.2, meaning that (8^2 — 8^1) ~ (1.2— l/v / 2)Ao, which is much larger than 
(0.754 — 1/ v / 2)Ao. If we go to much higher densities, where the calculation is under 
quantitative control, we find an even more striking enhancement: when g = 0.79 we find 
8^2/ Ao > 1000! We see that (relative to expectations based on experience with point- 
like interactions) the crystalline color superconductivity window is wider by more than 
four orders of magnitude at this weak coupling, and is about one order of magnitude 
wider at accessible densities if weak-coupling results are applied there. 2 

We have found that 8^2/^0 diverges in QCD as the weak-coupling, high-density limit 
is taken. Applying results valid at asymptotically high densities to those of interest in 
compact stars, namely /j ~ 400 MeV, we find that even here the crystalline color su- 
perconductivity window is an order of magnitude wider than that obtained previously 
upon approximating the interaction between quarks as point-like. The crystalline color 
superconductivity window in parameter space may therefore be much wider than previ- 
ously thought, making this phase a generic feature of the phase diagram for cold dense 
quark matter. The reason for this qualitative increase in 8/J2 can be traced back to the 
fact that gluon exchange at weaker and weaker coupling is more and more dominated 
by forward-scattering, while point-like interactions describe s-wave scattering. What is 
perhaps surprising is that even at quite large values of g, gluon exchange yields an order 
of magnitude increase in 8/J2 — 8//1 . 

This discovery has significant implications for the QCD phase diagram and may 
have significant implications for compact stars. At high enough baryon density the 
CFL phase in which all quarks pair to form a spatially uniform BCS condensate is 
favored. Suppose that as the density is lowered the nonzero strange quark mass induces 
the formation of some less symmetrically paired quark matter before the density is 
lowered so much that baryonic matter is obtained. In this less symmetric quark matter, 
some quarks may yet form a BCS condensate. Those which do not, however, will have 



LOFF condensates have also recently been considered in two other contexts. In QCD with \i u < 0, 
Hd > and /j„ = — hj, one has equal Fermi momenta for u antiquarks and d quarks, BCS pairing occurs, 
and consequently a (fid) condensate forms [28, 29]. If — /j u and /Jd differ, and if the difference lies in the 
appropriate range, a LOFF phase with a spatially varying (fid) condensate results [28, 29]. The result 
of Ref. [18] that the LOFF window is much wider than previously thought applies in this context also. 
Suitably isospin asymmetric nuclear matter may also admit LOFF pairing, as discussed recently in Ref. 
[30]. Here, the interaction is not forward-scattering dominated. 



differing Fermi momenta. These will form a crystalline color superconducting phase 
if the differences between their Fermi momenta lie within the appropriate window. In 
QCD, the interaction between quarks is forward- scattering dominated and the crystalline 
color superconductivity window is consequently wide open. This phase is therefore 
generic, occurring almost anywhere there are some quarks which cannot form BCS 
pairs. Evaluating the critical temperature T c above which the crystalline condensate 
melts requires solving the nonzero temperature gap equation obtained from (15) as done 
in Ref. [16] for the case of a point- like interaction. As in that case, we expect that all 
compact stars which are minutes old or older are much colder than T c . This suggests 
that wherever quark matter which is not in the CFL phase occurs within a compact star, 
crystalline color superconductivity is to be found. As we discuss in the next section, 
wherever crystalline color superconductivity is found rotational vortices may be pinned 
resulting in the generation of glitches as the star spins down. 

Solidifying the implications of our results requires further work in several directions. 
First, we must confirm that pushing Fermi surfaces apart via quark mass differences has 
the same effect as pushing them apart via a 8/j introduced by hand. Second, we must 
extend the analysis to the three flavor theory of interest. And, third, before evaluating 
the pinning force on a rotational vortex and making predictions for glitch phenomena, 
we need to understand which crystal structure is favored. 



GLITCHES IN COMPACT STARS 

We do not yet know whether compact stars feature quark matter cores. And, we do not 
yet know whether, if they contain quark matter, that quark matter is color-flavor locked, 
meaning that quarks of all colors and flavors participate in BCS pairing, or whether the 
BCS condensate leaves some quarks unpaired. The lesson we take from the toy model 
analysis is that because the interaction between quarks in QCD is dominated by forward 
scattering, rather than being an s-wave point-like interaction, the difference in Fermi 
momenta between the unpaired quarks need not fall within a narrow window in order 
for them to form a crystalline color superconductor. 

We wish now to ask whether the presence of a shell of crystalline color supercon- 
ducting quark matter in a compact star (between the hadronic "mantle" and the CFL 
"inner core") has observable consequences. A quantitative formulation of this question 
would allow one either to discover crystalline color superconductivity, or to rule out its 
presence. (The latter would imply either no quark matter at all, or a single CFL-nuclear 
interface [10].) 

Many pulsars have been observed to glitch. Glitches are sudden jumps in rotation fre- 
quency Q. which may be as large as AQ./Q. ~ 10~ 6 , but may also be several orders of 
magnitude smaller. The frequency of observed glitches is statistically consistent with the 
hypothesis that all radio pulsars experience glitches [31]. Glitches are thought to orig- 
inate from interactions between the rigid neutron star crust, typically somewhat more 
than a kilometer thick, and rotational vortices in a neutron superfluid. The inner kilome- 
ter of crust consists of a crystal lattice of nuclei immersed in a neutron superfluid [32]. 
Because the pulsar is spinning, the neutron superfluid (both within the inner crust and 



deeper inside the star) is threaded with a regular array of rotational vortices. As the 
pulsar's spin gradually slows, these vortices must gradually move outwards since the ro- 
tation frequency of a superfluid is proportional to the density of vortices. Deep within the 
star, the vortices are free to move outwards. In the crust, however, the vortices are pinned 
by their interaction with the nuclear lattice. Models [33] differ in important respects as 
to how the stress associated with pinned vortices is released in a glitch: for example, the 
vortices may break and rearrange the crust, or a cluster of vortices may suddenly over- 
come the pinning force and move macroscopically outward, with the sudden decrease in 
the angular momentum of the superfluid within the crust resulting in a sudden increase 
in angular momentum of the rigid crust itself and hence a glitch. All the models agree 
that the fundamental requirements are the presence of rotational vortices in a superfluid 
and the presence of a rigid structure which impedes the motion of vortices and which 
encompasses enough of the volume of the pulsar to contribute significantly to the total 
moment of inertia. 

Although it is premature to draw quantitative conclusions, it is interesting to speculate 
that some glitches may originate deep within a pulsar which features a quark matter 
core, in a region of that core which is in the crystalline color superconductor phase. A 
full three-flavor analysis is required, first of all in order to check whether the qualitative 
conclusions reached in the two-flavor analyses done to date persist, and second of all 
in order to determine whether the LOFF phase is a superfluid. If the only pairing is 
between u and d quarks, this 2SC phase is not a superfluid [4], whereas if the LOFF 
pairing involves the s quarks, as seems likely, a superfluid is obtained [1, 4]. Henceforth, 
we suppose that the LOFF phase is a superfluid, which means that if it occurs within a 
pulsar it will be threaded by an array of rotational vortices. It is reasonable to expect that 
these vortices will be pinned in a LOFF crystal, in which the diquark condensate varies 
periodically in space. The diquark condensate vanishes at the core of a rotational vortex, 
and for this reason the vortices will prefer to be located with their cores pinned to the 
nodes of the LOFF crystal. 

A real calculation of the pinning force experienced by a vortex in a crystalline color 
superconductor must await the determination of the crystal structure of the LOFF phase. 
We can, however, attempt an order of magnitude estimate along the same lines as that 
done by Anderson and Itoh [35] for neutron vortices in the inner crust of a neutron star. 
In that context, this estimate has since been made quantitative [36, 37, 33]. For one 
specific choice of parameters [15], the LOFF phase is favored over the normal phase 
by a free energy Floff ~ 5 x (10 MeV) 4 and the spacing between nodes in the LOFF 
crystal is b = 7l/(2|q|) ~ 9 fm. The thickness of a rotational vortex is given by the 
correlation length ^ ~ 1/A ~ 25 fm. The pinning energy is the difference between the 
energy of a section of vortex of length b which is centered on a node of the LOFF 
crystal vs. one which is centered on a maximum of the LOFF crystal. It is of order 
E p ~ ^loff^ 3 ~ 4 MeV. The resulting pinning force per unit length of vortex is of order 
f p ~ Ep/b 2 ~ (4 MeV) / (80 fm 2 ). A complete calculation will be challenging because 
b < ^, and is likely to yield an f p which is somewhat less than that we have obtained 
by dimensional analysis. Note that our estimate of f p is quite uncertain both because 
it is only based on dimensional analysis and because the values of A, b and Floff are 
uncertain. (We have a reasonable understanding of all the ratios A/Aq, 8/j/Aq, q/Ao 



and consequently bAo in the LOFF phase. It is of course the value of the BCS gap 
Ao which is uncertain.) It is premature to compare our crude result to the results of 
serious calculations of the pinning of crustal neutron vortices as in Refs. [36, 37, 33]. It 
is nevertheless remarkable that they prove to be similar: the pinning energy of neutron 
vortices in the inner crust is E p 1 — 3 MeV and the pinning force per unit length is 
f p « ( 1 - 3 MeV) / (200 - 400 fm 2 ) . 

A quantitative theory of glitches originating within quark matter in a LOFF phase 
must await further calculations, in particular a three flavor analysis and the determination 
of the crystal structure of the QCD LOFF phase. However, our rough estimate of the 
pinning force on rotational vortices in a LOFF region suggests that this force may 
be comparable to that on vortices in the inner crust of a conventional neutron star. 
Perhaps, therefore, glitches occurring in a region of crystalline color superconducting 
quark matter may yield similar phenomenology to that observed. Were this to happen, 
we can hope that a more detailed analysis would reveal distinctions among observed 
glitches, with some better understood as conventional glitches, originating in the inner 
crust, and others better understood as glitches originating deep within the star, in quark 
matter in the crystalline color superconductor phase. This is surely strong motivation for 
further investigation. 
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